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Abstract 

We use the corrections to the Newton-Einstein secular precessions 
of the longitudes of perihelia w of some planets (Mercury, Earth, Mars, 
Jupiter, Saturn) of the Solar System, phenomenologically estimated as 
solve-for parameters by the Russian astronomer E.V. Pitjeva in a global 
fit of almost one century of data with the EPM2004 ephemerides, in 
order to put on the test the expression for the perihelion precession 
induced by an uniform cosmological constant A in the framework of 
the Schwarzschild-de Sitter (or Kottler) space-time. We compare such 
an extra-rate to the estimated corrections to the planetary perihelion 
precessions by taking their ratio for different pairs of planets instead 
of using one perihelion at a time for each planet separately, as done so 
far in literature. The answer is negative, even by further re-scaling by 
a factor 10 (and even 100 for Saturn) the errors in the estimated extra- 
precessions of the perihelia released by Pitjeva. Our conclusions hold 
also for any other metric perturbation having the same dependence on 
the spatial coordinates, as those induced by other general relativistic 
cosmological scenarios and by many modified models of gravity. Cur- 
rently ongoing and planned interplanetary spacecraft-based missions 
should improve our knowledge of the planets' orbits allowing for more 
stringent constraints. 
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1 Introduction 

Introduced for the first time by Einstein (1917) to allow static homogeneous 
solutions to the Einstein's equations in the presence of matter, the cosmo- 



logical constant A, which turned out to be unnecessary after the discovery 
of the cosmic expansion by Hubble (1929), has been recently brought back 
mainly as the simplest way to accommodate, in the framework of general 
relativity, the vacuum energy needed to explain the observed acceleration 
of the universe (Riess et al., 1998; Perlmutter et al., 1999). For the rela- 
tion between the cosmological constant and the dark energy see (Peebles 
and Ratra, 2003). For a general overview of the cosmological constant see 
(Carroll, 2001) and references therein. Theoretical problems concerning the 
cosmological constant are reviewed in (Weinberg, 1989). 

Since, at present, there are no other independent signs of the existence 
of A apart from the cosmological acceleration itself, attempts were made in 
the more or less recent past to find evidence of it in phenomena occurring on 
local, astronomical scales with particular emphasis on the precession^ of the 
perihelia u! of the Solar System's inner planets (Islam, 1983; Cardona and 
Tejeiro, 1998; Wright, 1998; Kerr et al., 2003; Kraniotis and Whitehouse, 
2003; Dumin, 2005; lorio, 2006; Jetzer and Sereno, 2006; Kagramanova et 
al., 2006; Sereno and Jetzer, 2006, 2007; Adkins et al., 2007; Adkins and 
McDonneh, 2007). 

2 The perihelion precession induced by a uniform 
cosmological constant and the confrontation with 
the data 

Starting from the radial acceleration (Rindler, 2001) 

Aa = ^AcV, (1) 

where c is the speed of light, imparted by an uniform cosmological constant 
A in the framework of the spherically symmetric Schwarzschild vacuum solu- 
tion with a cosmological constant, i.e. the Schwarzschild-de Sitter (Stuchlk 
and Hledk, 1999) or Kottler (1918) space-time, Kerr et al. (2003) by using 
the Gauss equations for the variation of the Keplerian orbital elements (Roy, 
2005) worked out the secular, i.e. averaged over one orbital revolution, pre- 
cession of the pericenter of a test-body induced by the cosmological constant 

^For other local effects (gyroscope precession, mean motion change, geodetic precession, 
gravitational red-shift, deflection of light, gravitational time-delay, doppler tracking of 
spacecraft on escape trajectories) induced by A see, e.g., (Kagramanova et al., 2006; 
Sereno and Jetzer, 2006). 
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A finding 

HA = if— )^^^- (2) 



n 



In it n = GM/a^ is the Keplerian mean motion of the planet moving 
around a central body of mass M, G is the Newtonian constant of gravitation 
and a and e are the semimajor axis and the eccentricity, respectively, of the 
test body's orbit. 

Here we wish to offer an alternative derivation of eq. (2) based on the use 
of the Lagrange perturbative scheme (Roy, 2005). The Lagrange equation 
for the pericentre is 



dw 1 ^(Vpcrt) Vl-e^g(Vpcrt) 

dt na^Vl - e2 tanz di na^e de 

where i is the inclination angle to the equator of the central mass and (Vpcrt) 
is the perturbing potential Vpcrt averaged over one orbital revolution. For 
the Schwarzschild-de Sitter spacetime the cosmologically-induced additional 
potential is (Kerr et al., 2003) 

Va = - JacV2. (4) 

D 

By evaluating eq. (4) onto the unperturbed Keplerian ellipse defined by 

r = a{l — ecos E), (5) 

where E is the eccentric anomaly, and integrating over one orbital period 
Pb = 27r/n by means of 

dt = (^^^) dE, (6) 
/ (l-ecosS)^d-E = 7r(2 + 3e2) (7) 

JO 

yields 

(Va) = -^AcV(2 + 3e2). (8) 

By inserting eq. (8) into eq. (3) one obtains just eq. (2). 

Jetzer and Sereno (2006), Sereno and Jetzer (2007), Adkins et al. (2007) 
and Adkins and McDonnell (2007) obtained, in different frameworks, the 
same result of eq. (2). Note that (w)^ oc ^/aF{l^^(?), where, for a uni- 
form A, the proportionality factor is common to all the bodies orbiting a 
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given central mass. Moreover, eq. (2) was obtained by using the standard 
radial isotropic coordinate which is commonly used in the Solar System 
planetary data reduction process to produce the ephemerides (Estabrook, 
1971), so that eq. (2) can meaningfully be used for comparisons with the 
latest observational determinations of the non-Newtonian/Einsteinian secu- 
lar precessions of the longitude of the perihelia^ w (Pitjeva, 2005a). Indeed, 
they were estimated by contrasting, in a least square sense, almost one 
century of data of different kinds with the suite of dynamical force mod- 
els of the EPM2004 ephemerides (Pitjeva, 2005b) which included all the 
standard Newtonian and Einsteinian dynamics, apart from just any exotic 
effects as the ones by A on both the geodesic equations of motion and of the 
electromagnetic waves. Thus, such extra-precessions of perihelia, estimated 
independently of our goal, account, in principle, for any unmodelled force 
existing in nature. 

Since the cosmological accelerated expansion yields A ~ 10~^^ cm~^, 
Kerr et al. (2003) concluded that the precession of eq. (2) is too small to 
be measured in the Solar System. lorio (2006), Jetzer and Sereno (2006), 
Kagramanova et al. (2006), Sereno and Jetzer (2006), Sereno and Jetzer 
(2007) Adkins et al. (2007) and Adkins and McDonnell (2007) used eq. (2) 
and the extra-precessions of the inner planets of the Solar System estimated 
by Pitjeva (2005a) to put constraints on A. In particular, Jetzer and Sereno 

(2006) , after working out the effect of A on the pericentre of a general two- 
body system with arbitrary masses in the standard post-Newtonian gauge, 
used various binary pulsar systems and planets of the Solar System one at a 
time separately; Sereno and Jetzer (2006) discussed the possibilities offered 
by future interplanetary ranging, especially for Pluto; Sereno and Jetzer 

(2007) and Adkins et al. (2007) obtained eq. (2) as a particular case in the 
framework of cosmological models with a non-null acceleration. Adkins and 
McDonnell (2007) worked out the perihelion precessions under the action of 
arbitrary central forces obtaining eq. (2) as a particular case. In all such 
cases— and also in previous analyses performed when only data for Mercury 
existed'' (Islam, 1983; Cardona and Tejeiro, 1998; Wright, 1998; Kraniotis 

^The longitude of perihelion m is defined, for orbits nearly equatorial like the Solar 
System's ones, as the sum of the argument of perihelion uj and of the longitude of the 
ascending node f2; the latter one, sensitive to the out-of-plane disturbing forces, is not 
affected by the entirely radial A— induced extra acceleration. 

^In fact, also the non-Newtonian perihelion precession of the highly eccentric (e = 
0.826) orbit of the asteroid Icarus was investigated (Shapiro et al., 1971; Sitarski, 1992; 
Zhang, 1994; Shahid-Saless and Yeomans, 1994), but it was never used for putting con- 
straints on A. We will not use it in the present analysis. 
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and Whitehouse, 2003)— the general scheme followed was to assume the plan- 
ets, or the pulsar systems considered, separately one at a time and to derive 
constraints on A for each of them by considering it as a free parameter. 



3 Taking the ratio of the periheha 

Here we will follow, instead, a different approach which is able to tell us 
something much more definite about eq. (1) and eq. (2). We will construct 
the ratios of the estimated extra-precessions of perihelion for different pairs 
of planets A and B and will compare them to the corresponding ratios of the 
precessions of eq. (2) for the same planets. Thus, it is possible to construct 
the following quantity 



If eq. (2) is correct, independently of the value of A provided that, of course, 
it is nonzero and small enough to assure that the perturbation approach 
followed to derive eq. (2) is appropriate, then eq. (9) must be compatible 
with zero within the errors in w and a (Pitjeva, 2005b). It must be noted 
that such an approach holds, in general, also for any other extra- acceleration 
term, whatever origin it may have, of the form 



so that our conclusions will not be restricted to the Schwarzschild-de Sit- 
ter spacetime only. Indeed, there are other general relativistic cosmological 
models (Mashhoon et al., 2007) and many long-range models of modified 
gravity (Cognola et al., 2005) which are able to induce a perturbing accel- 
eration like that of eq. (10) (Ruggiero and lorio, 2007). 

3.1 The inner planets 

Let us start with the inner planets whose relevant orbital parameters are 
listed in Table 1. The answer they give is negative. Indeed, for the pairs 
A=Mars, B=Mercury and A=Earth, B=Mercury we have 




(9) 



(10) 



AMarMer = 7.8 ± 0.2, 



< 



(11) 



, AsarMer = 4.1 ± 0.2; 
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Table 1: Estimated semimajor axes a, in AU (1 AU= 1.49597870691 x 10^^ 
m) (Pitjeva, 2005b), and phenomenologically estimated corrections to the 
Newtonian-Einsteinian perihelion rates, in arcseconds per century (" cy~^), 
of Mercury, the Earth and Mars (Pitjeva, 2005a). Also the associated errors 
are quoted: they are in m for a (Pitjeva, 2005b) and in " cy~^ for w (Pitjeva, 
2005a). For the semimajor axes they are the formal, statistical ones, while 
for the perihelia they are realistic in the sense that they were obtained from 
comparison of many different solutions with different sets of parameters and 
observations (Pitjeva, private communication 2005). However, the results 
presented in the text do not change if 5uj are re-scaled by a factor 10. 



Planet 


a (AU) 


5a (m) 


w {" cy-^) 


Sw {" cy-i) 


Mercury 


0.38709893 


0.105 


-0.0036 


0.0050 


Earth 


1.00000011 


0.146 


-0.0002 


0.0004 


Mars 


1.52366231 


0.657 


0.0001 


0.0005 



a negative result at 40 — cr and 20 — 0" level, respectively. The other four 
pairs of inner planets yield, instead, results compatible with zero. The 
uncertainties in eq. (11) have been conservatively worked out by propagating 
the errors in'* w and a in eq. (9) and linearly summing the resulting biased 
terms 

the dominant sources of uncertainty are by far the perihelion rates. It is 
important to stress that even if bw were 10 times larger than the errors 
released by Pitjeva (2005a) and reproduced in Table 1— which are not the 
mere formal, statistical ones— we would still be forced to rule out eq. (2) at 
some a level. 

3.2 The outer planets 

It is interesting to use also some of the outer planets for which it was possi- 
ble to estimate the corrections to the perihelion precessions (Pitjeva, 2006), 
quoted in Table 2. The giant planets whose extra-precessions of the perihe- 
lion are at our disposal are Jupiter, Saturn and Uranus because the temporal 
extension of the data set used covered at least one full orbital revolution just 
for such planets: indeed, the orbital periods of Neptune and Pluto amount 

*The eccentricities are negligible. 
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Table 2: Estimated semimajor axes a, in AU (1 AU= 1.49597870691 x 10^^ 
m) (Pitjeva, 2005b), and phenomenologically estimated corrections to the 
Newtonian-Einsteinian perihelion rates (Pitjeva, 2006), in arcseconds per 
century (" cy~^), of Jupiter, Saturn and Uranus. Also the associated errors 
are quoted: they are in m for a (Pitjeva, 2005b) and in " cy~^ for w (Pitjeva, 
2006). For the semimajor axes they are the formal, statistical ones, while 
for the perihelia they are the formal ones re-scaled by a factor 10 in order 
to yield realistic estimates for them. 



Planet 


a (AU) 


6a (m) 


^ {" cy-') 


6w {" cy-i) 


Jupiter 


5.20336301 


639 


0.0062 


0.036 


Saturn 


9.53707032 


4,222 


-0.92 


2.9 


Uranus 


19.19126393 


38,484 


0.57 


13 



to about 164 and 248 years, respectively. For the external regions of the So- 
lar System only optical observations were used, apart from Jupiter (Pitjeva, 
2005b); they are, undoubtedly, of poorer accuracy with respect to those used 
for the inner planets which also benefit of radar-ranging measurements, but 
we will show that they are accurate enough for our purposes. Let us stress 
that in Table 2 we re-scaled by a factor 10 the otherwise formal, statisti- 
cal errors in the estimated extra-rates of perihelia. The pair A=Jupiter, 
B=Saturn yields 

AjupSat = 0.41 ± 0.06; (13) 

it is incompatible with zero at about 7 — a level. If we further re-scale by 
10 the error of Table 2 in the extra-precession of Saturn, for which Pitjeva 
did not use radiometric data, i.e. by 100 its formal, statistical error, we get 

AjupSat = 0.4±0.2, (14) 

negative at 2 — 0" level. In regard to Jupiter, since (Pitjeva, 2005b) used also 
some radiometric data for it, we believe that a re-scaling of 10 of the formal 
error in its estimated perihelion extra-rate is adequate. 
For the pair A=Jupiter, B=Mercury we have 

AjupMcr = 51 ± 12 (15) 

by using the figure of Table 1 for the uncertainty in the extra-precession of 
Mercury; a negative result at more than 4 — a level. If, in a very pessimistic 
approach, we re-scale by 10 the error of Table 1 in the Mercury's extra- 
rate— although it is not the formal one— we get 

AjupMcr = 51±34, (16) 
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which is incompatible with zero at 1.5 — a level. The errors in the semimajor 
axes, even if re-scaled by 10 or more, do not affect at all our results. The 
other pairs of planets yield results compatible with zero. 

4 Conclusions 

In this paper we used the latest observation-based determinations of the 
non-Newtonian/Einsteinian precessions of the longitudes of perihelia vj of 
some planets of the Solar System obtained by E.V. Pitjeva (Institute of Ap- 
plied Astronomy, Russian Academy of Sciences) (Pitjeva, 2005a, 2006) with 
the EPM ephemerides (Pitjeva, 2005b) to constrain the dynamical effects in- 
duced in our planetary arena by a non-zero uniform cosmological constant in 
the framework of the Schwarzschild-de Sitter (or Kottler) space-time. Such 
corrections to the standard perihelion precessions were determined without 
modelling at all the effects of A on both the geodesic equations of motion 
of planets and electromagnetic waves carrying information on them, so that 
they fully account, in principle, for such effects. The ratios of different 
pairs of planetary perihelia were used; by conservatively treating the errors 
in the estimated extra-precessions it turns out that the expression of the 
A— induced perihelion precession is ruled out at many a level. It is impor- 
tant to note that our phenomenological approach is quite general because 
it holds also for any small extra- acceleration of the form A = Ar, ^ : 
such a functional form is, in fact, induced not only by general relativistic 
cosmological models (Mashhoon et al., 2007) but also by modified models 
of gravity (see, e.g. (Cognola et al., 2005; Ruggiero and lorio, 2007)). 

The present analysis relies only upon the extra-precessions estimated 
so far by Pitjeva; it would be useful if also other teams of astronomers 
would estimate independently their own corrections to the standard plane- 
tary perihelion rates by exploiting the huge records of modern observations 
currently available. In regard to future perspectives, our knowledge of the 
motion of the inner planets of the Solar System should improve in the near 
future thanks to the ongoing Messenger (http:/ /messenger .jhuapl.edu/) and 
Venus Express (VEX) (Fienga et al., 2008) missions to Mercury and Venus, 
respectively, and to the planned BepiColombo (http://sci.esa.int/science- 
e/www/area/index.cfm?fareaid=30) (Milani et al., 2002) mission to Mer- 
cury; also Planetary Laser Ranging (PLR) (Chandler et al., 2005), e.g. 
to a lander on Mars (Merkowitz et al., 2007) or to the Mercury laser al- 
timeter (Sun et al., 2005), would greatly increase the accuracy in plane- 
tary orbit determination. In regard to the outer regions of the Solar Sys- 
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tern, maybe the processing of the radiometric data from the ongoing Sat- 
urnian mission of Cassini and from the Jupiter's flyby of New Horizons^ 
(http://pluto.jhuapl.edu/), recently occurred, might improve our knowledge 
of the motion of the outer planets as well in a not too far future. 
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